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Magnetic resonance imaging (MRI) is the cor¬ 
nerstone technique for diagnostic medicine, biol¬ 
ogy, and neuroscience [l]. This imaging method 
is highly innovative, noninvasive and its impact 
continues to grow [2]. It can be used for measur¬ 
ing changes in the brain after enhanced neural 
activity [ 3 ], detecting early cancerous cells in tis¬ 
sue | 4 ], as well as for imaging nanoscale biological 
structures [5], and controlling fluid dynamics [6], 
and it can be beneficial for cardiovascular imag¬ 
ing ra- The MRI performance is characterized by 
a signal-to-noise ratio, however the spatial reso¬ 
lution and image contrast depend strongly on the 
scanner design [5]. Here, we reveal how to ex¬ 
ploit effectively the unique properties of metasur¬ 
faces for the substantial improvement of MRI ef¬ 
ficiency. We employ a metasurface created by an 
array of wires placed inside the MRI scanner un¬ 
der an object, and demonstrate a giant enhance¬ 
ment of the magnetic field by means of subwave¬ 
length near-field manipulation with the metasur¬ 
face, thus strongly increasing the scanner sensi¬ 
tivity, signal-to-noise ratio, and image resolution. 
We demonstrate experimentally this effect for a 
commercially available MRI scanner and a biolog¬ 
ical tissue sample. Our results are corroborated 
by measured and simulated characteristics of the 
metasurface resonator, and our approach can en¬ 
hance dramatically functionalities of widely avail¬ 
able low-field MRI devices. 

It is well known that the strength of a signal deliv¬ 
ered by conventional MRI devices depends on the static 
magnetic field (Bo) of the imaging system [8j. Over 
the last two decades, low-field MRI machines operat¬ 
ing at Bo = 1.5 T have been used widely as key clini¬ 
cal tools. Recently, newly constructed high-field systems 
(Bo = 3 T) are being successfully exploited in hospi¬ 
tals ^ and there is a growing demand for ultra-high-field 
MRI machines (7 T or more) [TO] . 

A drive for high-field magnets is fueled by the benefits 
of higher signal-to-noise ratio (SNR) and better image 
resolution m- Still, to improve the overall image qual¬ 
ity one has to overcome the increasing static and radio¬ 
frequency (RF) magnetic field inhomogeneities that lead 
to the susceptibility artifacts and contrast variations jT2j . 
Further, the most important problem at high static fields 


5 0 > 7 T is a potential harm to a human body [13]. In 
particular, a possibility for tissue heating negligible at 
magnetic fields around 1.5 Tesla becomes substantial at 
higher fields due to an increase of the RF energy absorp¬ 
tion m- Also, an exposure to magnetic fields above 2 T 
can produce effects like vertigo and nausea, therefore the 
examinations above 3 T are conducted under especially 
careful medical supervision. Moreover, there exists quite 
a number of implants and medical devices that can be 
safe for low fields (1.5 T) and become unsafe at 3 T sys¬ 
tems US IE]. Therefore, the problem of enhancing the 
image quality without increasing the static magnetic field 
is very relevant to the MRI technology. 

Up to now, many different approaches have been pro¬ 
posed in order to enhance the MRI characteristics, while 
keeping the static field unchanged. The first way is the 
coil optimization. The progress in the RF coil technol¬ 
ogy over the last decades has already made an essential 
contribution on the MR scanner design. Comprehensive 
work has shown that it is possible to reduce scanning 
time with the help of parallel imaging methods [16, IT?] , 
that a substantial improvement in the SNR is available 
through multi-channel coils DEI, and that larger area can 
be examined using traveling-wave technique [H[20]. The 
second approach is based on the idea to employ spe¬ 
cial contrast agents [213 locally enhancing the magnetic 
field, such as rare-earth magnetic atoms [22] and mag¬ 
netic nanoparticles [23]. The third approach relies on 
special additional pads placed near the object. This can 
be done in several ways: by using high-permittivity pads 
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Figure 1. Artist’s view of a biological object placed on a 
metasurface resonator. The shaded rectangular region marks 
the region of interest. 
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in order to improve the RF field homogeneity and to re¬ 
duce the specific absorption ratio [24]; by using dielec¬ 
tric resonators (e.g. the cylindrical resonator when it is 
used as a coil with a hole in the center for an object) 
to increase the RF magnetic field and the SNR locally 
[25l 26 . However there exist several limitations for ap¬ 
plications of the dielectric resonators in MRI technology. 
These limitations include large geometrical dimensions of 
resonators at low frequencies, difficulties with placing an 
object at the maximum of the magnetic field inside the 
resonator, and difficulty to enhance only the RF magnetic 
field at the surface of the resonator, which is necessary 
for investigating large objects. 

The MRI efficiency can be further increased by using 
metamaterials, artificial structured media characterized 
by the values of effective permittivity and permeabil¬ 
ity not commonly available in Nature. Metamaterials 
demonstrate many fascinating fundamental properties, 
and they are very promising for engineering the prop¬ 
erties of electromagnetic devices and for tailoring the 
electromagnetic field-matter interactions [271 [28]. MRI 
is a very important field for applications of metamate¬ 
rials. Several existing implementations of metamaterials 
in MRI have been already proposed, such as Swiss-roll 
arrays [29], split-ring arrays with negative magnetic per¬ 
meability [30, 31 , magnetoinductive waveguides [32 , and 
endoscopes from arrays of metallic wires [33]. However, 
we believe that the full potential of the metamaterials for 
MRI technique is far from being uncovered. 

Here we demonstrate how 1.5 T MRI scanning effi¬ 
ciency is improved when the wire metamaterial-based 
surface (metasurface) structure is placed inside the scan¬ 
ner. Metasurfaces were studied by many groups [34] [35] 
and have shown spectacular applications in the long- 
wavelength regime for communication m and in optics 
to control the propagation of light and to provide sub¬ 
wavelength imaging [36l 137] . We propose a metasurface, 
based on the resonant array of metallic wires, that is 
inserted into the water based holder [Fig. [I]. We demon¬ 
strate that this metasurface can controllably redistribute 
the RF electromagnetic field at deep subwavelength level 
and simultaneously increases the SNR more than twice. 
The immediate result of this increase is that a higher res¬ 
olution image can be obtained over the same time slot or 
faster examination can be performed with similar reso¬ 
lution. Such metasurface presents a substantial step to¬ 
wards the improvement of MRI systems that are widely 
available in hospitals, without increasing the static mag¬ 
netic field Bq. 

Our main results and experimental setup are summa¬ 
rized in Fig. [2] We have realized a metasurface as an ar¬ 
ray of 14 x 2 wires. The studied biological sample (ex vivo 
fish) has been placed on the metasurface structure, em¬ 
bedded inside the water phantom, see the inset of Fig. IK 
and placed in a 1.5 T MR scanner, that has operating 
frequency 63.8 MHz [Fig. HR The birdcage coil, which 
resides in the bore of the scanner, has been used as a 
transmitting antenna see Supplementary Materials and 


Methods for more details. The radio frequency signal 
has been received by a surface body array. The scans 
have been performed according to the standard mag¬ 
netic resonance imaging protocol (see Methods). Pan¬ 
els (b)-(e) summarize the obtained results. Fig. It and 
Fig. If present the images obtained without the meta¬ 
surface structure for the long scanning time of 1020 s 
(b) and for the short time of 120 s (c). The SNR ratio 
is clearly increased for longer scanning time due to the 
accumulation of the signal. Panels (d) and (e) present 
the scans obtained for the fish placed on the metasur¬ 
face for the long (1020 s) and short (120 s) scanning time 
correspondingly. The amount of received valuable sig¬ 
nal has substantially increased with the metamaterial [cf. 
Fig. [ 2 ) 3,0 and Fig. [2]d,e] . Moreover the amount of received 
signal for the short time with metasurface is comparable 
with that for longer time without the metamaterial [cf. 
Fig. [ 2)3 and Fig. [2^]. Hence, Fig. 1 indicates, that it is 
possible to achieve higher signal and better image qual¬ 
ity at lower scan time. Due to the fact that MRI image 
quality is evaluated by SNR, we have performed more 
detailed experiments and quantitative numerical simula¬ 
tions, discussed below. 

First, we have made experimental measurements of 
the SNR in the region of interest (see the black plane 
of Fig. [T]) under the metasurface. Based on results of 
the separate scans of signal and noise in two configu¬ 
rations (phantom with the metasurface with and w/o 
the RF pulse and empty phantom with and w/o the 
RF pulse) we determined the ratio SNR 2 /SNRi « 2.7 , 
where SNR 2 corresponds to the ratio with the metasur¬ 
face inside the phantom and SNRi corresponds to the ra¬ 
tio for the empty phantom (see Supplementary Materials 
for more details). The increase of the SNR by more than 
a factor of two effectively corresponds using an MR sys¬ 
tem with a two times higher static magnetic field. Here, 
it is possible to improve MRI characteristics just putting 
the sample on the metasurface, without using the high 
field scanners. 

The enhancement of the signal with addition of meta¬ 
material (Fig. e) is due to the resonant coupling to 
electromagnetic modes of the metasurface. The length 
of each wire is tuned to the Fabry-Perot condition for 
the first eigenmode at the operating frequency of 1.5 T 
MRI machine, i.e. / = c/(^/eL) ~ 63.8 MHz. Here, 
81 is the dielectric constant of the background media, 
water in our case, and L = 25.5 cm is the wire length. 
For the first Fabry-Perot mode the highest magnetic field 
is localized in the middle part of the surface and the elec¬ 
tric field is localized near the edges of the wires [38j |39] . 
Hence, the observed effect can be explained by the RF 
magnetic field enhancement on the wire metasurface. In 
particular, the SNR ratio is determined as a ratio of the 
RF magnetic field in the scanned region (a signal) 
to the square root of the absorbed power (a noise) m 
The Fabry-Perot resonance leads not only to the overall 
enhancement of both B field and noise in the phantom, 
but it pins also the magnetic field to the scanned region 
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Figure 2. (a) A photo of the experimental setup placed into a MRI scanner. The inset illustrates schematically a metasurface 

embedded in the phantom and a sample. (b,c) MRI images (false colors) scanned without a metasurface for 1020 s and 120 s, 
respectively. (d,e) MRI images scanned with the metasurface for 1020 s and 120 s, respectively. 


suppressing the electric field and, hence, the absorbed 
power in this region. Such a relative redistribution of the 
spatial fields leads to a growth of SNR. It should be men¬ 
tioned that the metasurface enhances also the transmit¬ 
ted RF magnetic field B+ (see Supplementary Materials 
and Methods for more details). 

To further confirm that the MRI enhancement of the 
SNR is mediated by the Fabry-Perot modes of the meta¬ 
surface, we have compared experimentally and numer¬ 
ically two prototypes with different wire lengths L = 
25.5 cm and L = 20 cm. In order to determine the 
working frequency of metamaterial in numerical simu¬ 
lation we have placed a small magnetic loop as a source 
under the metasurface in the center and analyzed the 
source reflection coefficient Sn. The loop was located 
at 1.2 cm distance from the surface of the metamate¬ 
rial. The minimum of the reflection coefficient is reached 
at the metasurface resonance, when the power is effi¬ 
ciently transmitted to the metasurface. In our case this 
takes place at the frequency of the first eigenmode of 
the metasurface. Red and blue curves in Fig. [3] corre¬ 
spond to the reflection coefficient of the antenna for dif¬ 
ferent wire lengths, L = 25.5 cm (blue) and L = 20 cm 
(red). The Fabry-Perot mode frequency and the min¬ 
imum of the reflection coefficient exhibit red shift for 
longer wires. The Fabry-Perot resonance of the struc¬ 
ture with longer wires is tuned to the nuclear magnetic 
resonance at / = 63.8 MHz. This explains the observed 
enhancement of the MRI image, as shown in the middle 
column of Fig. [3 }d. For shorter wires, the Fabry-Perot 
resonance is detuned. The image quality measured for 
this structure is much weaker (see the right column in 
Fig. §>) and very similar to the results of the scanning 
of the fish without a metasurface (see the left column in 
Fig. [ 3 ) 3 ) . The subtle difference between the images with¬ 
out a metamaterial and with a short metasurface can 


be explained due to the excitation of different transverse 
modes of a finite metamaterial sample at the frequencies 
below the Fabry-Perot resonance m- The excitation 
of these modes leads to the redistribution of the elec¬ 
tromagnetic field and a slight enhancement of SNR (see 
the right bottom panel in Fig. [3 Jd). For even shorter 
wires this effect will be minimized further at the 1.5 T 
operational frequency and no difference between scans 
with short wires and without metamaterial is observed. 
Therefore, the detuned metasurface does not influence 
the imaging process. Thus, Fig. It conclusively demon¬ 
strates, that the experimentally observed MRI enhance¬ 
ment is due to the spectral matching between the first 
Fabry-Perot mode of the metasurface and the frequency 
of the nuclear magnetic resonance. 

The safety aspects in MRI are mainly determined by 
the specific absorption ratio, that quantifies the absorp¬ 
tion of electromagnetic energy and therefore the risk of 
tissue heating due to the application of the RF pulses, 
necessary to produce the MR signal. The specific ab¬ 
sorption ratio is proportional to the square of the induced 
electric field. The electric field on the surface of the meta¬ 
material is concentrated at the edges of the wires (see the 
absolute value of electric field induced on the top of the 
metasurface Fig. [4| and has an area of minimum values 
in the region of interest. Therefore, as soon as the object 
is placed on the resonator in a proper way the examina¬ 
tion is safe and even safer than without resonator. In 
our case the optimal region coincides with the region of 
interest (black rectangular plane in Fig. [l]) that is equal 
to 67% of the resonator length (see Supplementary Mate¬ 
rials for detailed discussion). This metasurface resonator 
provides the possibility to homogeneously enhance and 
redistribute the RF field in such a way, that the magnetic 
field is located in the area of the object and the electric 
field is removed from the region of interest. Moreover, 































4 



Figure 3. (a) Spectral dependence of the reflection coefficient 
for the magnetic loop antenna placed above the metasurface 
in the top plane center. Red and blue curves present the sim¬ 
ulated spectra and correspond to the two metamaterial struc¬ 
tures with different wire length, L — 20 cm and L = 25.5 cm, 
respectively, (b) Experimental MRI scans and numerically 
calculated SNR maps for three cases: without metamaterial, 
with metasurface length L = 25.5 cm and L — 20 cm respec¬ 
tively. 


due to the low electric field region the SNR enhancement 
does not depend on the different tissue permittivities of 
various bodies. Also, the enhancement principle differs 
qualitatively from the earlier demonstrated metamaterial 
designs, mainly used to transmit an RF signal from the 
MRI tube or as a matching device between a coil and 
a patient [nisniES]. Here, we enhance locally only the 
magnetic field in the scanned region while keeping the 
electric field and the SAR low. Although a metasurface 
made of split-ring resonators [30] could be employed po¬ 
tentially in the similar regime, it suffers from stronger 
electric field and more inhomogeneous magnetic field in 
the imaged area. Hence, the proposed structure is more 
advantageous for a number of reasons. 

With a view to underline the advantages of the meta¬ 
surface resonator we have compared the proposed struc¬ 
ture with the conventional high-permittivity resonator. 
The intrinsic property of the solutions to the Maxwell 
equations in the homogeneous medium is the coincidence 
between the antinodes of the magnetic field and the nodes 
of the electric one. Hence, it is impossible to maximize 
the RF magnetic field at the surface of the dielectric res¬ 
onator while keeping the electric field away from the ob¬ 
ject under examination. We have numerically studied 
a bulk rectangular dielectric resonator with parameters 
tuned to obtain the first mode at the same frequency 
/ = 63.8 MHz. In order to excite the resonator we have 
placed a magnetic loop antenna in the center near the 


surface (see bottom left corner of Fig. [4|. For qualita¬ 
tive comparison we have increased the number of layers 
in the metamaterial, so that the bulk dielectric resonator 
and the metamaterial resonator have approximately the 
same dimensions and are excited identically (see left cor¬ 
ner of Fig. |4|. Importantly, the metasurface behaves as 
a bulk wire metamaterial already for two or more layers 
of wires m ■ 

As can be seen in Fig. [4j for a wire metamaterial res¬ 
onator (top panels) the magnetic field is always localized 
at the resonator surface, while for a dielectric resonator 
(bottom panels) the magnetic field is localized inside the 
high permittivity structure, and electric field is localized 
at the surface. That is why the dielectric resonator can be 
employed only in the case when the object under study 
is placed inside of a high-permittivity resonator, which 
seems to be quite challenging. 

We have realized an unique metamaterial resonator in 
the form of a wire metasurface. The Fabry-Perot reso¬ 
nance of the wire array is spectrally tuned to the nuclear 
magnetic resonance. This allows us to observe a sub¬ 
stantial enhancement of the signal-to-noise ratio enabling 
a reduction of the MRI scanning time and also obtain 
higher image resolutions. Also it provides an opportunity 
to locally control the electric and magnetic fields near 
the object. There exists quite a number of possible im¬ 
plementations of metasurfaces for medical examinations 
inspired by the concepts of smart clothes and wearable 
electronics (see Supplementary Materials for details). For 
example, metal strips can be printed on clothes, so that a 
patient will have to wear a special lightweight jacket be¬ 
fore the examination. In this particular case, the meta¬ 
surface produces a homogeneous excitation from all sides. 
This involves no health risks provided that the strip ends 
are isolated from the patient, while it gives a possibility 
to use any desirable radio-frequency coils positioned near 
the patient. The metasurface can be embedded in cer¬ 
tain parts of the patient table as well. To summarize, the 
use of metasurface is a promising pathway to improve the 
MRI efficiency and design simple yet effective scanning 
protocols. 


METHODS 

All experiments were performed using a 1.5 T whole- 
body MRI system (Philips Healthcare) in combination 
with the 16 channel Torso XL receive array. The fish 
was located inside the water phantom and the receive 
array was placed around the water phantom (Fig. [2^i) . 
Two identical measurements (with the same position of 
the receiving coil) were performed: (1) a fish inside the 
water phantom without a metasurface; (2) a fish inside 
water phantom on a metasurface. 

The metasurface was made from nonmagnetic brass 
wires. Wire length L, radius r and lattice period a are 
equal to 25.5 cm, 0.1 cm and 1 cm, respectively. The 
wires were supported by a double-layer plastic film, fill- 
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Figure 4. Comparison between electric and magnetic field distributions in two resonators: wire metamaterial (on the top) and 
dielectric resonator (on the bottom). Left row illustrates the high-permittivity and metamaterial resonators that have been 
used in numerical simulations. Top line shows the schematic cut planes. 


able by water. The period of the structure was optimized 
to achieve a homogeneous magnetic field enhancement in 
the region of interest. The number of layers was chosen 
to reach the highest magnetic field enhancement, while 
keeping the thickness of the structure as thin as possible. 

Initially, we have determined by the direct measure¬ 
ments of the B 0 map that there is no influence of the 
metasurface on the main static magnetic field of MRI. 
The scans shown in Fig. [2] were obtained using a stan¬ 
dard spin-echo sequence, with the following parameters: 
field of view was 300 mm by 300 mm, flip angle was 
90°, the echo time (TE) was 20 ms, for long time scans 
(Fig. |2]3,d) the repetition time (TR) was 1240 ms and 
for short scan time (Fig. [2j^,e) the TR was 480 ms. For 
SNR measurements we obtained separate data for signal 
and noise. The measurement sequence was carried out 
twice: the first time with the RF pulses and therefore 
acquiring signals and noise and the second time without 
any RF pulses and therefore acquiring only noise. The 
signal was defined as the mean pixel intensity value in a 
region of interest, while the noise was defined as the stan¬ 
dard deviation in this region in the noise image. For this 
measurement we used gradient-echo sequence (GRE), the 


flip angle was 15°, the TR was 50 ms and TE was 3.7 
ms. After an addition of the metamaterial, the RF power 
levels that previously generated 90 ° excitation pulses in 
the region of interest created a much larger angle due 
to enhancement of the RF field. In order to avoid over¬ 
tipping with metasurface present, the power level of the 
excitation pulses were re-optimized (reduced). 

Numerical simulations were performed using the time- 
domain solver of CST Microwave Studio 2013 package. 
The small brass ring, excited by the discrete port (see 
red triangle in the left schematic plots of Fig. [4| was 
used as a source. 
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